Conductance characteristics between a normal metal and a two-dimensional 
Fulde-Ferrell-Larkin-Ovchinnikov superconductor: the Fulde-Ferrell state 
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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state has received renewed interest recently due 
to the experimental indication of its presence in CeCoIns, a quasi 2-dimensional (2D) d-wave super- 
conductor. However direct evidence of the spatial variation of the superconducting order parameter, 
which is the hallmark of the FFLO state, does not yet exist. In this work we explore the possibility 
of detecting the phase structure of the order parameter directly using conductance spectroscopy 
through micro-constrictions, which probes the phase sensitive surface Andreev bound states of d- 
wave superconductors. We employ the Blonder- Tinkham-Klapwijk formalism to calculate the con- 
ductance characteristics between a normal metal (N) and a 2D s- or d x 2_ y 2-wave superconductor 
in the Fulde-Ferrell state, for all barrier parameter z from the point contact limit (z = 0) to the 
tunneling limit (z> 1). We find that the zero-bias conductance peak due to these surface Andreev 
bound states observed in the uniform d-wave superconductor is split and shifted in the Fulde-Ferrell 
state. We also clarify what weighted bulk density of states is measured by the conductance in the 
limit of large z. 

PACS numbers: 74.25.Fy, 74.50.+r 



C 

O 

Q : 

(N 
> 

o 

in 
o 

-I— > 



o 
o 



X 



I. INTRODUCTION 



In the early 1960's, Fulde and Ferrelli and Larkin and 
Ovchinniko-^ proposed the possibility that a supercon- 
ducting state with a periodic spatial variation of the or- 
der parameter would become stable when a singlet su- 
perconductor is subject to a large Zeeman splitting. The 
Zeeman splitting could be due to either a strong mag- 
netic field or an internal exchange field. Under such a 
strong magnetic or exchange field, there is a splitting of 
the Fermi surfaces of spin-up and -down electrons, and 
condensed pairs of electrons with opposite spins across 
the Fermi surface may be formed to lower the free en- 
ergy from that of a normal spin-polarized state. These 
pairs have a non-zero total momentum 2q, which causes 
the phase of the superconducting order parameter to vary 
spatially with the wave number 2q. This state is known 
as the Fulde-Ferrell (FF) state. Larkin and Ovchinnikov 
(LO), on the other hand, proposed independently an al- 
ternative scenario, in which the order parameter is real, 
but varies periodically in space, possibly in more than one 
directions. Both types of states are now known (collec- 
tively) as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state. It has not yet been observed in conventional 
low-T c s-wave superconductors, which are mostly three- 
dimensional, probably because the orbital effect of the 
magnetic field dominates the Zeeman effect. The situ- 
ation has been changed by experimental results sugges- 
tive of the FFLO state in heavy-fermion, quasi- ID or- 
ganic, or high-T c superconductorsj&liSiSiL&iLiSiAi many 
of these compounds are quasi-one or two-dimensional, 
thus the orbital effect is weak when the magnetic field is 
in the conducting plane or along the chain. Recent ex- 



perimental results in CeCoIns, a quasi-2D d-wave super- 
conductor, are particularly encouragingiiSiiiii^iiSiiSiiLiSi 

We note in passing that this subject is also of interest to 
the nuclear and particle physics communities because of 
the possible realization of the FFLO state in high den- 
sity quark matter and nuclear matter^ as well as in cold 
fermionic atom systems^ 

The most unambiguous evidence of the FFLO states 
should be based on phase-sensitive experiments that can 
directly reveal the spatial variation of the phase or sign 
of the order parameter. One possibility is the Josephson 
effect i2i To the best of our knowledge there has been no 
report on this or other phase sensitive measurements thus 
far. In this paper, we consider an alternative phase sen- 
sitive probe, i.e., the conductance spectroscopy through 
a micro-constriction. A powerful method to calculate the 
differential conductance [G(V) = dI(V)/dV] character- 
istics of a normal metal/superconductor junction (NSJ) 
was developed by Blonder, Tinkham, and Klapwijk 22 , 
which unified quasiparticle tunneling and Andreev refl.ee- 
tiort 2 ^ 2 ^. Only s-wave superconductors were considered 
in that work, but the method has since been general- 
ized to the ci-wave casei 25 : 26 When the sign of the or- 
der parameter experienced by the electron- and hole-like 
quasiparticles (QPs) before and after specular reflection 
at the junction interface changes, zero-energy Andreev 
bound states (ABSs) (also called the midgap states^) 
are formed at the S side of the N/S interface (for z ^ 0). 
These ABSs will give rise to a zero-bias conductance peak 
(ZBCP) in the tunneling spectrum of the NSJ. A well 
known example of such NSJ is when the junction inter- 
face along the nodal line of the d x 2_ y 2-wave superconduc- 
tor [(110) contact] i 25 i 26 i 27 i 28 i 29 i 30 This feature has since 



2 



been widely used to identify the order-parameter symme- 
try of unconventional superconductors. 31 ' 32 ' 33 - 34 - 35 ' 36,37 
We would like to emphasize here that the ABSs are conse- 
quences of the phase change of the <i-wave order parame- 
teo 3 ^ 3 ^, and thus their spectra should also be sensitive to 
the spatial variation of the order parameter. In this pa- 
per, we will show this sensitivity by explicit calculations, 
and from their spectra detected via conductance spec- 
troscopy through a micro-constriction, one can extract 
the momentum of the superconducting order parameter. 
In the present work, we will focus on the FF state; the 
more general LO cases are currently being investigated 
and will be presented elsewhere As we will show be- 
low, for a d x 2_ y 2-wa,ve superconductor in the FF state, 
the ZBCP observed in a (110) contact is split and shifted 
by both the Zeeman field and pair momentum; the latter 
can be determined from the splitting. 

This paper is organized as follows. In section II we in- 
troduce the model and present the self-consistent mean- 
field solutions for both s- and ci :c 2_ y 2-wave superconduc- 
tors in the FF phases. The electron density of states 
(DOS) of the FF states is also obtained. The scheme 
to calculate the conductance characteristics is presented 
in section III, with numerical results for both s-wave 
and d x 2_ y 2-w&ve cases, and their relation with the cor- 
responding bulk electron DOS is discussed. The features 
associated with the ABSs in the FF phase are explored 
in section IV, where we solve their spectra in the FF 
state. We present further discussion on our results and 
a summary in section V. In this work, we consider zero 
temperature only, and calculations on d-wave supercon- 
ductors are simply referred to d x 2_ y 2-wave. 



II. THE FULDE-FERRELL STATE AND BULK 
DENSITY OF STATES 

We begin with the mean field Hamiltonian for a 2D 
superconductor in the BCS state (q = 0) or the FF state 
(q^O): 

~Hmf = + cr A i o-B)cj io .Ck CT 

kcr 

- 5Z( Ak q C k+qT C -k+ql + A kqC-k+qlC k+qT ), (1) 

k 

where £ k is the single-particle kinetic energy relative to 
the Fermi energy ep, B is the Zeeman magnetic field, fj,Q 
is the magnetic moment of the electron, A kq is the pair- 
ing potential and satisfies the self-consistent condition, 



A k q - - Vkk'(C-k'+qXCk'+qT) 



(2) 



For s-wave, we have Vkk' = — Vo, while for d-wave, we 
have Vkk' = -V cos(20 k ) cos(20 k -) (here, |£ k |, |f k /| < 
huj c <C £f, w c is the cutoff, 9k is the azimuthal angle of 
k). The orbital coupling between the magnetic field and 



electron motion is absent when the field is in the plane. 
The mean field Hamiltonian could be rewritten as 



k 



k+ql 



const. (3) 



where 



£k+q + —Akq 



-A 



kq 



-6 



k-q 



/i -B 



To diagonalize it, we perform the Bogoliubov-Valatin 
transformation 



C k+qT 
C -k+ql 



u k Vk W a k 



(4) 



and choose 



«k _ £ k s) + £k 



u k 



^kq 



As) 
2k_ 

Ek 



l-\Vk\\ 



where Ek = yj A£ q + ^ 2 , = (£ k+q + £k_ q )/2, and 

£k^ = (£k+ q — 6c- q )/2, from which we get the diagonal- 
ized Hamiltonian 

Hmf = ^(£ k +a k "k + # k _/3 k Ac) + const. (5) 

k 



with eigenenergies (ct = ±1) 

Ek* = Ek + <r( t i B + (,l a) ). 



(6) 



There are regions in the k-space where the Cooper pairs 
are destroyed and occupied by electrons of one spin 
species; these are states with Eka < 0. In these cases 
the Bogoliubov-Valatin transformation (Eq. Q) should 
be replaced by 



C k+qT 
-k+ql 

when Ek+ < 0, or 



-v£ Uk 



Ck+qT \ _ ( " k <'k 



(7) 



(8) 



'-k+ql / \^"k «k / \ Pk 

when Ek- < 04i Then, the diagonalized Hamiltonian is 
Hmf = ^2(\Ek+\alak + |£fc_|j8£&) + const. (9) 

k 

with positive quasiparticle energies. 

In the weak coupling limit, the total energy of the sys- 
tem is given by 

f 2^ s) \vk\ 2 , £k,±i>0l A 2 
(H - /iV> = < Ck+q + MoS, Ek+ < > - TT- 
k I Ck-q-MoS, E k - < I 



(10) 
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FIG. 1: Qualitative sketch of the phase diagram of an s-wave 
superconductor under a Zeeman field, including the possibil- 
ity of Fulde-Ferrell state. Solid line shows the phase transition 
between the superconducting state and the normal spin po- 
larized state; dash-dotted line indicates the phase transition 
between the usual BCS state and the FF state; dotted line 
gives the Clogston-Chandrasekhar critical field. 



Here, for an s-wave superconductor, we have Ak q = 
A q , and for a d-wave superconductor, we have Ak q — 
Ag cos 20k- For a given Zeeman field B, we calculate the 
energy of the pairing state numerically using Eq. I|1(J|) 
with Q = qvF /Aq varied, and compare it with that of the 
normal spin polarized state, in order to find the ground 
state. In our numerical calculation presented below, we 
take hoj c — 10A (here A is the gap of the usual BCS 
state at T = without a Zeeman field, and for <i-wave it 
is the gap along antinodal direction or maximum gap), 
and H = hqB/Aq. The calculation of the self-consistent 
pairing potential shows that the superconducting state 
is destroyed at H > 1 for s-wave^ and at H > 1.06 for 
d-wavem^^ In the superconducting regime, for s-wave 
case, at H p i ps 0.704, a transition occurs and the FF 
state becomes the ground state. For d-wave case, the 
critical Zeeman field is H p i ps 0.544 where an FF phase 
occurs with q along the nodal direction. At a higher 
H' pl w 0.78, the FF state with q along the antinodal di- 
rection dominates^ The critical fields H p i are slightly 
smaller than the Clogston-Chandrasekhar fields which is 
H p i = l/v2 for s-wave^, and « 0.56 for d-wave^ (see 
also Ref. I46T) . The qualitative illustrations of the phase 
diagrams are presented in Figs. Q and |21 

The finite momentum of the order parameter leads to 
a nonzero supercurrent in the ground state. However, 
the supercurrent is expected to be compensated by the 
backflow of the quasiparticle current, and therefore, we 
have a zero total current i Here we demonstrate that the 
total current is indeed zero for the self-consistent mean- 
field ground state. The total current operator can be 
written as 

J = ~ J2(( k + q)4+q T C k+qT - ( k - q) C -k+qi C -k+ql) 

k 

(11) 
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FIG. 2: Qualitative sketch of the phase diagram of a d-wave 
superconductor under a Zeeman field, including the possibil- 
ity of Fulde-Ferrell state. The curves have the same meaning 
as in Fig.0 except that the dashed line is the phase boundary 
separating FF states with q along the antinodal direction and 
the nodal direction. Here, a is the angle between q and the 
antinodal direction. 

and in the superconducting state, its expectation value 
is 

f 2qM 2 , E K±1 > 1 

(J) = J2{ q + k, £ k+ <0 }. (12) 

m k [ q - k, £ k < . J 

By differentiate Eq. I|10(l at fixed A we can verify that 

(J)=-e^-(H-tiN). 

Since (Ti. — /iN) is minimized in the q, A space, the cur- 
rent must be zero. 

The electron DOS can be evaluated using 

P {E)= Yj \u*\ 2 {5{E-E^ + ) + 8{E-E^)) 

k,B k± >0 

+ (\^\ 2 S(E + E k+ ) + \u k \ 2 5(E~E^)) 

k,_E k+ <0 

+ (\^\ 2 S(E-E k+ ) + \v k \ 2 5(E + E^)). 

k,_B k _<0 

(13) 

Within the approximation of ~ £k arid ^ k a ' ) w 
qvF cos(9\ L — 9q), it can be seen that for an arbitrary state 
k, there is always another state k' with £k' = — £k> so 
that both states have the same energy, and their weight- 
ing factors in Eq. 1131 add up to unity. Thus, Eq. 113f) can 
be simplified to 

p[E) = PMj\ e dC(6(E-\E k+ \) 
+ 5{E-\E k -\)), 
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FIG. 3: Representative electron density of states for various 
Zeeman fields and both s- and d-wave superconductors, (a) 
s-wave: solid, H = 0.704, Q = 0.804; dash, H = 0.8, Q = 
0.892; dot, H = 0.9, Q = 0.952. (b) d-wave with q along the 
antinodal direction: solid, H — 0.78, Q = 0.808; dash, H = 
0.88, Q — 0.9. (c) d-wave with q along the nodal direction: 
solid, H = 0.544, Q = 0.608; dash, H = 0.68, Q = 0.816; dot, 
H = 0.776, Q = 0.94. 



where p n (0) is the DOS of the normal state at the Fermi 
level. Finally, we obtain 



where 



Pe(E, 



Pe{E) _ 1 

p n (0) 2ir J 



e\ \A(k F ) 



d0p e (E,( 



(14) 



|A(k F ) 



with e a = E — a(poB + qvp cos(#k — ^q))- The result is 
presented in Fig. 



mean field Hamiltonian Q), we obtain 

Eu a (-x) = (h + ap B)u a (-x.) 

+ J dx'A(s,r)u CT (x'), (15a) 
Ev a (x.) = -(ho + <7p B)v a (x.) 

+ J rfx'A*(s,r)u CT (x'), (15b) 

where s = x — x', r=(x + x')/2, h — — V 2 /2m + 
US(x) - p and A(s,r) = / dk e ik s A(k, r)e i2c > r . Here, 
q is determined by minimizing the ground-state energy 
[Eq. (|10|) ]. To avoid accommodating super/normal cur- 
rent conversion at the N/S interface, which may require a 
nontrivial modification of the order parameter structure 
near the interface, we assume that q is parallel to the N/S 
interface at x = 0; this choice makes theoretical analysis 
simpler and comparison to experiments easier. Neglect- 
ing the proximity effect at the N/S interface, we have 
A(k, r) = A(k)0(x), where 0(a;) is the step function, 
and A(k) is the pairing potential in the bulk supercon- 
ductor. In the WKB approximation, the BdG equations 
have the special solutions of the form, 



( u <y\ = ,k r .x/ e lc i x u CT \ 



(16) 



where u,y and v a obey the Andreev equations 



eu a = „ *( k f + q) . VM g + A(k F )9(x)f} CT , (17a) 



m 

i(k F - q) 



• Vu ff + A*(k F )6(a;){t CT , (17b) 



where e = E — crpoB — q • kp jm. The eigenenergy E is 
symmetric about E = crp Q B + q • kp/m instead of zero 
as in Ref. These equations are similar to those of 
Ref. 0j where the authors studied conductance charac- 
teristics in the presence of a supercurrent along the junc- 
tion (note that the q/A° in Ref. is equal to Q/2 of 
this work) . Following the method of Blonder-Tinkham- 
Klapwijk 22 * 2 ^^, the normalized tunneling conductance 
at zero temperature is given by G = G ns /G nn , where G ns 
and G nn are the tunneling conductances in the supercon- 
ducting and normal states, respectively. Since the QPs 
of the two spin species are uncoupled, the conductance 
will simply be the average over the two spin components: 



G), 



1 



G„ 



III. THE CONDUCTANCE CHARACTERISTICS 

In the presence of a normal metal/superconductor 
interface, to study the QPs, one needs to solve the 
Bogoliubov-de Gennes (BdG) equations,— From the 



and the same for G nn . For simplicity, we will drop the 
spin index when there is no ambiguity. The conductances 
of each spin species are given by 



e 



tt/2 



G ns = / d6(l + \a(-E, 

n J-tt/2 



\b(E,e)\ 2 ), (is) 



5 



tt/2 



tt/2 



d6(i- \b{+oo,e)\ 2 ), 



where 



a(E,( 



cos" 



b(E,( 



?7 + (cos 2 9 + z 2 ) — rj- 
z(z + i cos 6)(i] + — 



V-) 



77+ (cos 2 9 + z 2 ) — rj_z 2 



e±^/e 2 -\A(9, 



= E — (JHqB — 



qk Fy 



(19) 



(20a) 



(20b) 



(20c) 



(20d) 



Here, a(E,9) and b(E,9) are the Andreev and normal 
reflection coefficients, respectively; 9± = 9 ± a, 9 is the 
angle between kp and the +a; axis, a is the angle between 
the antinodal direction and the +x axis for d-wave, and 
zero for s-wave, z = mU /kp is a dimensionless barrier- 
strength parameter. We give the normalized tunneling 
conductance of the FF states in both s- and c?-wave su- 
perconductors with the arrangement of q parallel to the 
N/S interface in Figs.EJEL an d also show the spin split- 
ting effect of the tunneling conductance at a large z = 5.0 
in Fig. [7| In Fig. [HJ we give the conductance of the com- 
peting uniform BCS states at the critical fields for both 
pair types of superconductors. Note that the Zeeman 
splitting for the s-wave case shown in Fig. |SJa) repro- 
duces the well-known results reviewed in Ref. |4£| and the 
z = results given in Ref. IHoL 

We start our discussion with an s-wave superconduc- 
tor or a d-wave superconductor with (100) contact, where 
there are no ABSs; thus the conductances are determined 
by the bulk quasiparticles. We note that, in the arrange- 
ment considered here, due to the effect of pair momen- 
tum, the conductance curves at large z seem to no longer 
coincide with the corresponding electron DOS for either 
case. [Compare Figs 0Jc),|5fc), and Etc) with Figs. El (a), 
(b), and (c).] At first sight, this appears to be against 
our intuitive understanding of what the NS J conductance 
at high-barrier limit (tunneling limit) is supposed to be 
measuring; we now resolve this issue below. 

In the arrangement of NSJs considered here, the ap- 
plied voltage and the measured current are both along 
a fixed direction in the conducting plane that is normal 
to the N/S interface. Thus we realize that the junction 
conductance in the high-barrier-strength limit is actu- 
ally measuring a cos 2 0-weighted DOS. That is, the QPs 
of various momenta k on the 2D Fermi surface (i.e., cir- 
cle) do not all make equal contributions to the junction 
conductance, but should be weighted by cos 2 9 where 9 
is the angle between k and the current direction x. The 
weighted DOS [p w (E)] measured in the high-z-limit junc- 
tion conductance is therefore: 



p w {E) 

Pn(0) 



1 



d9 p e (E, 9)\ cos 6f 



(21) 



In the uniform BCS state (without involving ABSs), such 
weighted average simply returns to the original DOS, i.e., 




FIG. 4: The normalized conductance vs voltage for normal- 
metal/s-wave FF superconductor junction: (a) z — 0, (b) 
z = 0.5, (c) 2 = 5.0. Solid, H = 0.704, Q = 0.804; dash, 
H = 0.8, Q = 0.892; dot, H = 0.9, Q = 0.952. 
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FIG. 5: Same as Fig.^Jexcept it is now for normal-metal/d- 
wave FF superconductor junction with (100) contact. Solid, 
H = 0.78, Q = 0.808; dash, H = 0.88, Q = 0.9. 
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FIG. 6: Same as Fig.^Jexcept it is now for normal-metal/d- 
wave FF superconductor junction with (110) contact. Solid, 
H = 0.544, Q = 0.608; dash, H = 0.68, Q = 0.816; dot, 
H = 0.776, Q = 0.94. 



FIG. 7: The normalized conductance G (solid line) and 
its two spin components G (dashed line for spin-up, dash- 
dotted line for spin-down) at a large z — 5.0. (a) s-wave 
with H = 0.704; (b) d-wave with q along antinodal direction 
and H = 0.78; (c) d-wave with q along nodal direction and 
H = 0.544. 



Eq. (|21|l is the same as Eq. (|14fl in the cases studied here. 
This is because the order parameter is isotropic over the 
momentum space for an s-wave superconductor, while for 
an N/(d-wave superconductor) junction with (100) con- 
tact, the order parameter is symmetric about 8 — ±7r/4 
so that the partial DOS of a fixed direction is also sym- 
metric about the same line. Because the cos 2 6 weighting 
factor adds up to be 1 for two angles that are symmetric 
about this line, we thus have the weighted DOS by high- 
z junction conductance to be again the same as the un- 
weighted DOS< 26 ^ However, the situation is very different 
in the FF state, because this kind of symmetry is broken 
by the pair momentum which causes a ^-dependent en- 
ergy shift, thus the two kinds of DOS are no longer the 
same. This is illustrated in Fig. [5] It is clear from this 
figure that the high-z junction conductance measures the 
cos 2 ^-weighted DOS, and not the un-weighted DOS in 
general. (The slight discrepancy between the tunneling 
conductance and the weighted DOS is because z — 5.0 is 
still not high enough.) 

The most prominent features of the high-z junction 
conductance of an N/(ci-wave FF superconductor) with 
(110) contact are due to the ABSs, which are the main 
focus of the present work. To interpret these features 
we need to understand how the pair momentum affects 
the spectra of the ABSs, which is the subject of the next 
section. 
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FIG. 8: The normalized conductance vs voltage for the com- 
peting uniform BCS states at the critical fields with differ- 
ent barrier intensities z. (a) s-wave superconductor with 
H — 0.704; (b) d-wave superconductor with H = 0.544 and 
(110) contact. 



IV. ANDREEV BOUND STATES IN THE 
D-WAVE FULDE-FERRELL 
SUPERCONDUCTOR WITH (110) JUNCTION 

As has been reviewed in the introduction, for an N/(d- 
wave superconductor) junction with (110) contact, a 
ZBCP is expected due to the formation of the "zero en- 
ergy Andreev bound states" (or "midgap states" ) at the 
junction interface. Thus for an N/(d-wave superconduc- 
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FIG. 9: Comparison of the high-z(=5.0) junction conduc- 
tance Gns (circles with a thin dotted line through them) with 
the un-weighted (dashed line) and cos 2 ^-weighted (solid line) 
DOS for (a) an N/(d-wave BCS superconductor) junction 
with (100) contact; (b) an N/(s-wave FF superconductor) 
junction, at H = 704; and (c) an N/(d-wave FF supercon- 
ductor) junction with (100) contact and q along (010), at 
H = 0.78. 



tor in an FF state) junction with (110) contact and q 
along (110), we need to understand the effects of [i B 
and q on the spectra of ABSs before we can understand 
the conductance characteristics for this junction. 

In the limit of z — * oo we can consider a simple model 
with an FF superconductor in the region x > 0, and vac- 
uum or an insulator in the region x < 0, with q parallel 
to the N/S interface. In the superconductor region, from 
Eq. (|17(l , we can find the solutions of bound states (which 
decays to zero as x — ► +oo) to be of the form, 



( " CT± ^ = e~ lx ( " CT± 1 

V V a ± ) ~ \ Va± J ' 



where 



y = mV\&\ 2 -e 2 /\k Fx \, 



u„±/va± = A/(e T iy/\A\ 2 -e 2 ) , 

and e = E — o~HqB — qkp y /m (kp x — kpeosO, kp y = 
kp sin 9) . In order to satisfy the vanishing wave- function 
boundary condition at x — 0, we need to first superpose 



the above solutions of opposite k x : 



J <y+ 



-■"iv i 



imposing the vanishing boundary condition at x = then 
yields 



Vn 



Mo 
Vn 



(22) 



For an s-wave superconductor, and for a d-wave super- 
conductor with (100) contact, Eq. I|22l) yields no ABSs 
solutions. For a d-wave superconductor with (110) con- 
tact, we obtain one solution for each — kp < kp y < kp 
at e = or 



E = ohqB + qkpy/m. 



(23) 



We see that the energies of the ABSs are first split by the 
Zeeman energy to ufJ-aB and then, shifted by an amount 
proportional to both the pair momentum q and the sine 
of the incident angle (which ranges between —90° and 
+90°); the combined results of them will lead to a split- 
ting, shifting, and broadening of the ZBCP, which we 
now analyze. With (110) contact, the order parameter 
is proportional to sin 29, which vanishes at 9 = 0, and 
±90°, implying that near these special values there are 
either no ABSs, or their contributions to conductance 
will be very weak because these are very loosely bound 
states; the dominant contributions come from ABSs with 
9 around ±7r/4 (or around gap maxima); the energy shift 
of these states due to pair momenta are in opposite direc- 
tion when the sign of 9 are different. Thus for sufficiently 
large q, which is the case here, the junction conductance 
at positive bias should exhibit two peaks, one on each side 
of /iqB, and a dip at the Zeeman field energy E = fioB. 
(The conductance has a symmetry about zero bias in the 
approximation adopted here, so we do not need to con- 
sider negative bias.) For the two peaks, the one on the 
right side (E p+ ) arises from < 9 < ir/2, whereas the 
one on the left side (E p -) arises from —ir/2 < 9 < 0. If 
only spin-up QPs are considered, the two peaks are of 
equal strength and equal distance from the dip as illus- 
trated in the dashed line of Fig. EJc). When the contri- 
butions from the QPs of both spin species arc summed 
up, the two peaks will not be symmetric about the dip 
(E p - will be shifted slightly to the right, but the influ- 
ence on E p+ and the dip is negligible), and a weak peak 
at zero bias emerges (see Fig. Etc)). To locate the po- 
sition of the peak E p+ , we first notice from Eas. ll8tf2Tj| 
that the bias voltage difference between the peak and the 
dip, 5E p = \E p+ — fioB\, should be a function of qv F in 
the high-z limit and vanishes when the pair momentum is 
zero. Numerically, therefore we can consider a simplified 
situation where there exists the pair momentum without 
Zeeman field (such in Ref. l48j) and calculate this differ- 
ence as a function of the pair momentum. The result is 
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15.0 




FIG. 10: The bias voltage difference of the tunneling con- 
ductance peak (E = E p +) and the dip (E = fioB) at large 
z(=20), SE P = \E P + — noB\, as a function of qvF (both in 
unit of A q ) for d-wave superconductor with (110) contact. 
The value range of qu^/Ag is extended to 0-20, while the 
physical range is about 1.5-7.7. Solid line is the linear fit of 
the numerical data, and its slope is 0.667 ± 0.001. The in- 
set is an illustration of the voltage difference measured in the 
experiment. 

shown in Fig. By fitting the data to a straight line 
through the origin, we obtain 

2 

E p+ « li B + -qv F . (24) 

Thus, by measuring the bias voltages of the peak and the 
dip in the high-z junction conductance with (110) con- 
tact, and in particular the difference between them, we 
can obtain a good estimate of q. We note that without 
this pair momentum, we would have all ABSs at energies 
E = <7fi B, which would have given rise to one sharp 
peak only at eV — fi^B in the same conductance plot, as 
shown in Fig. Efb) . Therefore we conclude that the sig- 
nature of the FF state is clearly revealed in the junction 
conductance characteristics, especially at high z, but the 
conductance behaviors at low z for all three cases stud- 
ied here are also quite novel, since they are quite different 
from the corresponding results for uniform s- or d-wave 
superconductors. 

V. SUMMARY AND DISCUSSION 

In this paper we have studied the conductance char- 
acteristics of a junction between a normal metal and a 
superconductor in the Fulde-Ferrell (FF) state, using the 
Blonder-Tinkham-Klapwijk formalism. We have studied 
both s- and d-wave cases, and for the latter case, we 
considered junctions along both the nodal [(110)] and 
antinodal [(100)] directions. 

The conductance characteristics of a micro- 
constriction is presumably easiest to understand in 
the high tunneling-barrier limit, when the conductance 



should give information about DOS of the supercon- 
ductor. In the FF phase, the Zeeman field should split 
the contributions to the conductance by the spin-up 
and -down QPs, i.e., shifting their contributions by the 
Zeeman energy in opposite directions. In addition, for 
QPs of either spin species, their contributions should be 
shifted by an amount proportional to the pair momen- 
tum, with a proportionality constant depending on the 
cosine of the angle between its kinetic momentum and 
the pair momentum. This proportionality constant can 
range from a negative maximum to a positive maximum. 
Thus we found the effect of q to be a broadening 
rather than a shift. One would expect similar effects 
to occur on the contributions by the ABSs, resulting 
in the shift and broadening of the ZBCP. Thus, in 
the high-barrier- limit, one might (naively) expect the 
tunneling spectrum of d-wave superconductor with (110) 
contact to be composed of broadened ZBCP's centered 
around ±fj,oB. In our numerical results, we find instead 
that the high-barrier-limit tunneling conductance of 
d-wave superconductors with (110) contact has a dip 
at the Zeeman energy, with one round peak on each 
side of it, and also another weak peak at exactly zero 
energy. This is quite different from the situation of 
the BCS state in the presence of a Zeeman field, in 
which case one expects the sharp ZBCP to be shifted to 
±/iqB, as well as the naive expectation above. The dip 
and the two round peaks can be understood as due to 
the fact that the q values appearing here are so large 
that they are already beyond the critical value obtained 
in Ref. \4§L which studied directly the current effect 
on the conductance characteristics in the absence of a 
Zeeman field. For such large q values, their effect on the 
ZBCP in the high-barrier-limit tunneling conductance 
is to split the ZBCP into two round peaks with a 
center dip. It is worth noting that such high values 
of q are not accessible through direct application of 
a supercurrent. The weak peak at zero energy turns 
out to be the result of summing up spin-up and -down 
contributions. Furthermore, numerical analysis shows 
that the bias voltage difference between the dip and the 
round peak on its right side is proportional to the pair 
momentum and thus gives us a simple way to estimate 
the pair momentum for d-wave superconductor with 
(110) contact. 

For s-wave superconductor and d-wave superconduc- 
tor with (100) contact, there is no ABS and the con- 
ductance is due to contributions from bulk quasiparti- 
cles exclusively. In these cases we have also found con- 
ductance features in the FF superconductors that are 
very different from the BCS superconductors. We found 
that because of the energy shift due to the pair momen- 
tum (which breaks the spatial symmetries in the original 
system), the conductance in the high barrier limit is no 
longer the same as the electron DOS; instead, it reflects 
a directionally-weighted DOS. In principle, by compar- 
ing the conductance and the bulk DOS 46 i 51 i 52 i 53 that are 
measured by other means (such as tunneling along the c- 
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direction instead of an in-plane direction discussed here) , 
one can also distinguish between BCS and FF states. 

The FF state studied here is the simplest version of the 
general Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) super- 
conductors. Its simplicity allows for a fairly straight- 
forward calculation of the conductance characteristics 
through micro-constriction, 40 « as well as an understanding 
of the results. Study of the conductance characteristics 
for general FFLO states is currently under way and will 
be reported elsewhere. 
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